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1 Introduction 

Buoyancy effects in unstably stratified mixing layers express themselves through 
gravity currents of heavy fluid which propagate in an ambient lighter fluid. These 
currents are encountered in numerous geophysical flows, industrial safety and 
environmental protection issues Q . During transition to turbulence a strong 
distortion of the separating interface between regions containing 'heavy' or 'light' 
fluid arises. The complexity of this interface will be used to monitor the progress 
of the mixing. We concentrate on the enhancement of surface-area and 'surface- 
wrinkling' of the separating interface as a result of gravity-effects. 

Mixing in a turbulent flow may be visualized by a number of snapshots, e.g., 
in an animation. However, for a quantitative analysis, more specific measures 
are needed. Here we represent the separating interfaces as level-sets of a scalar 
field c (0 < c < 1). We establish a close correspondence between the growth of 
surface-area and wrinkling, in particular at large buoyancy effects. We also show 
that this process can be simulated quite accurately using large-eddy simulation 
with dynamic subgrid modeling. However, the subgrid resolution, defined as the 
ratio between filter-width A and grid-spacing h, should be sufficiently high to 
avoid contamination due to spatial discretization error effects. 

2 Surface-area and wrinkling of level-sets 

To specify the intricate structure of the separating interface we consider global 
properties of the level-set defined by S(t) = {x | c(x, t) — 1/2}. The mixing- 
efficiency rj is defined in terms of the surface-area A as: 

„(*) = ;A(t)=f dS = [ dx*(c(x,t)-i)|Vc(x,i)| (1) 

Js(t) Jv 1 

Instead of integrating directly over the complex and possibly fragmented inter- 
face S(t), we evaluate the volume integral infjL in which S C V, using a locally 
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exact quadrature method developed in Q . A further characterization of the 
geometry is obtained by considering the average curvature C and wrinkling W: 



C{t)= f dS K (x,t) ; W(t)= I dS\K{x,t)\ 
Js(t) Js(t) 



?(*) J s(t) 
where the local curvature n is defined as 



(2) 



«(x,t) = V.n ; n(x) = ^^ T (3) 

|Vc(x)| 

The effect of gravity is represented by source terms in the momentum and 
energy equations 0] . Within the Boussinesq approximation we have: 

dtp + dj(pUj) = 

d t {pui) + djiputiij) + dip - dj(Tij = -(^r^Jfe 

/ ape \ 

dte + dj((e + p)Uj)-dj(aijUi) + djqj = -\-^) u 2 (4) 

where p denotes the 'light' fluid density, Ui is the velocity in the Xi direction, p 
the pressure, e the total energy density, qj the heat-flux and <xy = Sij/Re with 
Sij the rate of strain tensor and Re the Reynolds number. The source terms 
on the right hand side represent gravity effects in the negative vertical (i.e. X2) 
direction in which a — {ph — p)/p with p^ the density of the 'heavy' fluid. The 
Froude number Fr = U* j ' \fg*L* with g* the gravitational constant and U* , L* 
a reference velocity and length scale respectively. We consider a — 1. Finally, c 
denotes a scalar field which relates to the local density p as p = (1 — c)p + cp%. 
The scalar field c is assumed to evolve according to 



dt(pc) + dj{u 3 pc) - Y c d " c = ( 5 ) 



where Sc denotes the Schmidt number. The scalar field c is convected by the 
flow field Uj but in turn affects the flow through the source terms in |. The 
initial state for c is a step from a value in the lower half to 1 in the upper half 
of the domain, representing a heavy fluid layer on top of a lighter one. 



3 Buoyancy effects in turbulent mixing 

To quantify buoyancy effects in turbulent mixing we consider the temporal mix- 
ing layer as documented in Q . The flow starts from a perturbed laminar profile. 
A fourth order accurate finite volume discretization was used. Rapid transition 
to turbulence emerges, displaying helical pairings. At large buoyancy effects the 
separating interface becomes very distorted and protruding heavy/light fluid 
'fingers' with a complex geometrical structure emerge, as shown in figure 0(a). 
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Figure 1: Separating interface at t = 40 using Fr = 2 (a); surface-area, wrinkling 
and curvature (b) at Fr = 2: A(t) - A(0) (solid), W (dashed), C (dash-dotted) 
and Fr = 5: A(t) - A(0) (*), W (o), C (+). Resolution of 128 3 cells; = 10. 



The geometrical properties A, C and W of the separating interface are com- 
pared in figure 0(b) for two different Froude numbers. The surface-area and 
wrinkling display approximately the same behavior, while the average surface- 
curvature remains quite limited, indicating that regions with negative and posi- 
tive local curvature are about equally important in this flow. 




Figure 2: Mixing-efficiency rj (a) and surface wrinkling (b) at Fr = 2 and Sc = 5 
(solid), 10 (dashed), 20 (dash-dotted), 50 (dotted), 100 (o). Resolution: 128 3 . 

In figure |2| w e compare 77 and W for various Sc. An increase in Sc reduces 
dissipation in |5] and results in enhanced mixing. After a laminar development 
up to t « 10, a nearly algebraic growth r\ ~ t" emerges with (3 = (3(Sc). The 
wrinkling also develops algebraically. Initially this development is independent 
of Sc. Beyond f w 10 the algebraic growth-rate increases with Sc. 

We also considered LES of this mixing flow. In figure p| we compare pre- 
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Figure 3: Mixing-efficiency r\ at Fr — 2 and Sc = 10: DNS (solid), no-model 32 3 
(dashed), dynamic mixed model (32 3 : dash-dotted, 64 3 : dotted); A = L/16. 

dictions from DNS and LES with the dynamic mixed model. The filter-width 
A = L/16 with L the side of the flow domain. The mixing-efficiency is slightly 
under-predicted on 32 3 cells. The discrepancy with DNS is largely of numerical 
nature. By increasing the subgrid-resolution r = A/h from 2 on 32 3 to 4 on 
64 3 a strong improvement arises even though A was kept fixed. The remaining 
errors are more representative of subgrid modeling deficiencies Q . 

4 Concluding remarks 

Buoyancy can strongly increase turbulent mixing. A reasonable agreement be- 
tween surface-area and wrinkling was observed. This may be relevant, e.g., in 
combustion modeling to capture flame-speed and intensity. Prediction of mixing- 
efficiency using LES with dynamic models, showed a large influence of spatial 
discretization errors. A subgrid-resolution r = 4 is suggested. 
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